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Abstract. Various theoretical and algorithmic aspects of inverse problems in discrete to- 
\Q • mography of planar Penrose model sets are discussed. These are motivated by the demand 

of materials science for the reconstruction of quasicrystalline structures from a small number 
. of images produced by quantitative high resolution transmission electron microscopy. 
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1. Introduction 

Discrete tomography is concerned with the inverse problem of retrieving information about 
some discrete object from (generally noisy) information about its incidences with certain query 
sets. A typical example is the reconstruction of a finite planar point set from its line sums 
| in a small number of directions. More precisely, for a direction u £ S 1 (the unit circle), the 

(discrete parallel) X-ray X U F of a finite subset F of the Euclidean plane IR 2 in direction u 
gives the number of points of the set on each line in IR 2 parallel to u, i.e., X U F is the function 
\Q ' X U F :C U — ► N := N U {0}, defined by 

O ' 

^ X U F(£) :=card(Fn^) = ^l F (x), 

| where C u is the set of lines in direction « in E 2 , with obvious generalization to higher di- 

mensions. In the classical setting, motivated by crystals, the positions to be determined form 
a subset of a translate of the square lattice 7Z? or, more generally, of arbitrary lattices L in 
H d , where d > 2. The cases d = 2 and d = 3 are practically relevant. In fact, many of the 
problems in discrete tomography have been studied on 1? , the 'classical planar setting' of 
discrete tomography; see p]|2l[HlE]. 

In the longer run, by also having other structures than perfect crystals in mind, one has 
to take into account wider classes of sets, or at least significant deviations from the lattice 
structure. As an intermediate step between periodic and random (or amorphous) Delone sets 
(defined below), we consider systems of aperiodic order, more precisely, of so-called model sets 
(or mathematical quasicrystals) , which are commonly accepted as a reasonable mathematical 
model for quasicrystalline structures in nature |E|. 

The main motivation for our interest in the discrete tomography of model sets comes 
from the question how to reconstruct three-dimensional (quasi) crystals, or planar layers of 
them, from their images under quantitative high resolution transmission electron microscopy 
(HRTEM) in a small number of directions. In fact, in ®E|, a technique called QUANTITEM 
(quantitative analysis of the information coming from transmission electron microscopy) is 
described, based on HRTEM, which can effectively measure the number of atoms lying on 
lines parallel to certain directions. In particular, with the growing importance of surface sci- 
ence, there is also increasing interest in additional methods for the reconstruction of planar 
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structures which can supplement STM approaches. At present, the measurement of the num- 
ber of atoms lying on a line can only be achieved for some crystals; see (HUH- However, it is 
reasonable to expect that future developments in technology will improve this situation. It 
seems thus timely to investigate the mathematical foundations now. 

Here, we restrict ourselves to an example, namely to the well-known class of planar model 
sets Ap C IR 2 that are associated with the well-known Penrose tiling, and present some 
results on the discrete tomography of these sets, with an emphasis both on reconstruction 
and uniqueness problems. Note that proofs are omitted; details and extensions will appear 
in (ElE]. 



2. Penrose Model Sets 

We always let Cs = e 27ri//5 , a primitive 5th root of unity in C Then, Q(Cs) is the correspond- 
ing cyclotomic field, an algebraic number field of degree 4 over Q, and ^[£5] is its subring of 
cyclotomic integers. 

Remark 1. Let C5 denote the cyclic group of order 5, i.e., C5 = Moreover, C5 is 

understood to be supplied with the discrete topology. It is well-known that every z € ^[£5] 
can uniquely be written as z = Y^j=o a j{ z )Cli with cij(z) 6 7L. Let 02 be the automorphism 
of the Galois extension Q(Cs)/Q that is given by Cs 1 — > Cf- Identifying IR 2 and C in the 
canonical way, 02 gives rise to a map 

.~ : Z[C 5 ] -^R 2 x (R 2 x C 5 ), 



defined by z > — > yz, (c2(z),J2j =0 (aj(z) (mod 5))) J. Via projection on the second factor, this 

induces a map * : Z[C 5 ] — > R 2 x C 5 , defined by z 1 — ► (a 2 {z), E^oK'^) (mod 5))). Then, 
Z[Cs]~ is a lattice in IR 2 x (IR 2 x C5), i.e., a co-compact discrete subgroup. In fact, is 
the Z-span of the set {1 , (£5) , (C|) , (Cf) }• Finally, note that Z[( 5 ]* is dense in IR 2 x C 5 ; 
see 

It is well known by now that model sets arise from so-called cut and project schemes, 
compare d21E3- m particular, the class of Penrose model sets (PMS) arises from the following 
cut and project scheme; cf. |14j : 

(1) 

R 2 < — IR 2 x (IR 2 x C 5 ) — ► IR 2 x C 5 

U dense U lattice U dense 
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Z[C 5 ] ^ {(z = ai(z),{a2{z),^(a j (z)(mod5)))^\z€Z[C 5 ]} ^ Z[C 5 ]* 

j=0 




Given any subset W C R 2 x C5 with / W° C W C W° and W° compact, a so-called 
window, and any t G R 2 , we obtain a planar model set A(t,W) := t + /1(W) relative to the 
above cut and project scheme £Q) by setting /1(W) := {z £ Z[Cs] | G W}. 

Let P be the convex hull of the 5th roots of unity, which is a regular pentagon centred at 
the origin. Set := P, W<® := -P, := tP and := -rP, with r = (1 + v^)/2 
the golden ratio, and 

4 

:= (J (w& x {j (mod 5)}) C R 2 xC 5 . 
3=1 

Moreover, for it G R 2 , set := (u,0 (mod 5)) + W P , (W$)® := (u,0 (mod 5)) + (W® x 
{j (mod 5)}) and 4g := A(W%). If 4£ is generic, i.e., if one has W£ D Z[C 5 ]* = 0, then all 
translates of yip, meaning the sets t + Ap with t G R 2 , are called Penrose model sets. Note 
that this formulation avoids the usual ambiguities from non-minimal embeddings into 5-space. 

Remark 2. A p is not generic, while generic examples are obtained by shifting the window, 
i.e., Ap is generic for almost all u G R 2 . Joining any two points with distance 1 in a generic Ap 
by edges results in a Penrose tiling, which is a tiling with two types of rhombi. See Figure Q 
for a generic example; different generic choices of u result in locally indistinguishable (LI) 
Penrose tilings. Note that Penrose model sets /lp C R 2 are aperiodic, meaning that they 
have no translational symmetries. Further, Penrose model sets are Delone sets, i.e., they are 
uniformly discrete and relatively dense; cf. |12j . 
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3. Problems in Discrete Tomography of Penrose Model Sets 

Let Ap be a PMS, k G N and R > 0. A finite subset C of Ap is called a convex set in Ap 
when its convex hull contains no new points of Ap, i.e., when one has C = conv(C) n Ap. We 
denote by F(Ap), F<k(Ap), V < n{Ap) and C(Ap) the set of finite subsets of Ap, the set of 
finite subsets of Ap having cardinality < k, the set of subsets of Ap with diameter less than 
R and the set of convex subsets of Ap, respectively. 

Remark 3. The uniform discreteness of Penrose model sets Ap immediately implies the 
inclusion V < ^(Ap) C ^{Ap). 

Clearly, in order to obtain electron microscopy images of high resolution, one should allow 
only directions which yield dense lines in Penrose model sets. These directions are clearly 
contained in the set of all directions, called Ap -directions, which are parallel to a non-zero 
element of the difference set of Ap, i.e., 

Ap-Ap := {A- A' | A, A' G A P } C Z[C 6 ] ■ 

Calling a direction u G S 1 a Z[Cs]- direction when it is parallel to an element of Z[Cs] \ {0}, one 
has the following result. 

Proposition 1. If Ap is a PMS, the set of Ap- directions equals the set of ^[C,^]- directions. 

Let us indicate the main algorithmic problems in discrete tomography of Penrose model 
sets. For a direction u G S 1 , we use £«^ 5 ' to denote the set of elements l u of C u that pass 
through a point of ^[Cb]- Let ui, . . . , u m G S 1 be m > 2 pairwise non-parallel Z [(^-directions. 
The corresponding consistency, reconstruction and uniqueness problems are defined as follows. 

Consistency. 

Given functions p Ui : C Ui — > No, i G {1, . . . ,m}, whose supports are finite 
and satisfy supp(p u J C Cul^, decide whether there is a finite set F which is 
contained in a PMS and satisfies X Ui F = p Ui , for all i £ {1, . . . , m}. 
Reconstruction. 

Given functions p Ui : C Ui — ► No, i G {1,... , m}, whose supports are finite 
and satisfy supp(p Ui ) C C^f^; in the case that Consistency is satisfied, 
construct a finite set F which is contained in a PMS and satisfies X Ui F = p Ui , 
for alH G {1, . . . , m}. 

Uniqueness. 

Given a finite subset F of a PMS, decide whether there is a different finite 
set F' that is also a subset of a PMS and satisfies X Ui F = X Ui F', for all 
i G {1,... ,m}. 

In general, the above problem Reconstruction can have many solutions of rather different 
shape. Therefore, one is also interested in uniqueness results, e.g., the (unique) determination 
of the set 

U 

A P PMS 

or suitable subsets thereof by the X-rays in a small number of Z[Cs]-directions. More precisely, 
we define the concept of determination and the interactive concept of successive determination 
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as follows. Let £ be a collection of finite subsets of R 2 and let U C S 1 be a finite set 
of directions. We say that £ is determined by the X-rays in the directions of U if, for all 
F, F' G £ , one has 

(X U F = X U F', VueU) F = F' . 

We say that £ is successively determined by the X-rays in the directions of £/ = {u\, ■ ■ ■ , u m }, 
if, for a given F G £ , these can be chosen inductively (i.e., the choice of Uj depends on all 
X Uk F with k G {1, . . . ,j — 1}) such that, for all F' G £, one has 

(X U F' = X U F, Vti G i7) =► F' = F . 

We say that £ is determined (resp., successively determined) by m X-rays if there is a set U of 
m pairwise non-parallel directions such that £ is determined (resp., successively determined) 
by the X-rays in the directions of U. 

4. Computational Complexity and Uniqueness Results 

Let us begin with a result on computational complexity, where we apply the real RAM- 
model of computation, see [Ej. Here, each of the standard elementary operations on reals 
counts only with unit cost. This leads to the following tractability result. 

Theorem 1. When restricted to two S[Cs]- directions, the problems Consistency, Recon- 
struction and UNIQUENESS can be solved in polynomial time in the real RAM-model. 

Remark 4. It seems to be rather obvious from the results in P] that one cannot expect a 
generalization of Theorem ^ to the case of three or more Z[Cs]-directions. More precisely, we 
expect that, when restricted to three or more Z[Cs]-directions, the problems Consistency, 
Reconstruction and Uniqueness are NP-hard. 

Let us now present results dealing with the (successive) determination of finite subsets of 
Penrose model sets. Though we are not interested in non-Z[(^5]-directions for practical reasons, 
we begin with the following observation. 

Proposition 2. If Ap is a PMS and u G S 1 is a non-Z^s]- direction, the class of finite subsets 
F{Ap) is determined by the single X-ray in direction u. 

This last result immediately follows from the observation that, for all PMS yip, each line in 
the plane in a non-Z[Cs]-direction passes through at most one point of Ap, the latter being the 
reason for the practical irrelevance of this result. On the other hand, the next result shows 
that any fixed finite number of X-rays in Z[Cs]-directions does not suffice to determine the 
whole class of finite subsets of a fixed PMS Ap. 

Proposition 3. Let Ap be a PMS and U C S 1 an arbitrary, but fixed finite set of pairwise 
non-parallel ^[£5]- directions. Then, the set J-(Ap) is not determined by the X-rays in the 
directions of U. 

In order to obtain results on uniqueness, one has to restrict the class of finite sets under 
consideration. Within the class of finite subsets of a fixed PMS Ap with bounded cardinality, 
there is the following result. 



6 



M. BAAKE AND C. HUCK 



Proposition 4. Let Ap be a PMS and k G N. T/ien, i/ie se£ T<k{Ap) is determined by 
any set of k + 1 pairwise non-parallel ^[£5]- directions, while any set of 1 + |l°g2 pairwise 
non-parallel X-rays in directions is insufficient for this purpose. 

This last result is once again of limited relevance in practice, because typical atomic struc- 
tures to be determined comprise about 10 6 to 10 9 atoms, and, in order not to damage or even 
destroy the examined structures, one has to make sure that one uses at most 4 or 5 X-rays. 

Proposition 5. Let Ap be a PMS and R > 0. Then, the set V < n(Ap) is determined by two 
X-rays in Z[Cs]- directions. 

Though the last result seems to be more satisfactory, it is probably still of restricted use 
in practice. Here, the reason is that, in general, the second Z[£s]-direction cannot be chosen 
in such a way that it yields dense lines in Penrose model sets Ap, in other words, one would 
have to deal with images of poor resolution. A deeper result is the following, which deals with 
the class of convex subsets of a fixed PMS yip. 

Theorem 2. There is a set U C S 1 of four pairwise non-parallel Z[Cs]- directions such that, 
for all PMS Ap, the set C(Ap) is determined by the X-rays in the directions of U, while, for 
all PMS Ap and any set U C S 1 of three or less pairwise non-parallel Z[Cs]- directions, the set 
C{Ap) is not determined by the X-rays in the directions of U . 

For example, the set of Z^J-directions parallel to the elements of the following set U has 
the desired property to determine C(/lp) by the X-rays in its directions, 

(2) U :={(! + t) + ( 5 ,(t- I) + ( 5 ,-t + ( 5 ,2t- ( 5 } • 

Remark 5. By a result of Pleasants these directions can yield dense lines in Penrose 
model sets. It follows that, in the practice of quantitative HRTEM, the resolution coming 
from these directions is rather high, which makes Theorem look promising. 

Above, we restricted the class of finite subsets of a fixed Penrose model set Ap under 
consideration. In order to obtain positive uniqueness results, a second option is to consider 
the interactive technique of successive determination. One has the following positive results. 

Theorem 3. // Ap is a PMS, the set ^(Ap) is successively determined by two X-rays in 
Z[Cs]- directions, while the set 

U ^p) 

Ap PMS 

is successively determined by three X-rays in Zj[(s]- directions. 

Unfortunately, this result is again somewhat limited in practice because, in general, one 
cannot make sure that all the Z[(^5]-directions which are used match dense lines in Penrose 
model sets. 

Final Remark 

For further details in this spirit, we refer to [21 01 EJ^] for the lattice case and ^3^] for 
cyclotomic model sets which also provide a systematic generalization of the setting explained 
here for the PMS. 
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